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Abstract 

The extent of plastic flow in a spherical solid (assumed to 
be homogeneous and elastically and plastically isotropic) surroundin 
a concentric rigid sphere as a function of applied external pressure 
is calculated. The applied pressure necessary to cause plastic 
deformation throughout the solid is obtained. 




INTRODU'-VnON 


When a solid matrix surrounding a rigid sphere is pressurized, 
the rigid particle acts as a stress concentrator. The radial press- 
ure -at the interface at the rigid sphere exceeds the external 
hydrostatic pressure. When the applied pressure is increased 
'Sufficiently, yielding of the matrix begins at the interface. This 
is a well known phenomenon and can be observed around inclusions 
in solids where dislocations are punched out . ^ The pressure which 
causes the onset of yielding is available from earlier calculation 
from linear elasticity theory and Tresca or von Mises yield 
criteria. In the present paper we calculate the extent of 
plastic flow versus applied pressure and finally obtain the 
applied pressure necessary to cause plastic deformation through- 
out the solid. The present analysis assumes that there is no 
strain hardening. 

In the present paper external pressurization of a medium 

surrounding a rigid spherical inclusion is studied. Two methods 

are used. The first, developed here, leads more easily to the 

2 

solution than the second, which is the technique used by Hill 
for studying the internal pressurization of a hollow sphere. 





THEORETICAL 

Consider the matrix of external radius b surrounding a rigid, 
sphere of radius a. The boundary conditions are that the radial 
component of stress is given by 


a = -p at r=b 
rr ext 


( 1 ) 


and the radial component of displacement is given by 


u = 0 at r = a 


C2) 


Elastic Behavior 

At sufficiently small values of the material behaves 

everywhere in a linear elastic manner. The stresses then liave 
the form 


a - A + —3 

rr r 


A 

2r^ 


(3) 

(4) 


and the radial displacement is given. by 




(S) 


Here is the circumferential strain, E is Young's modulus 

u y 

V is Poisson's ratio and we have used the fact that n, ,=o‘ 

(pq) 0 0 



From these five equations , the state of stress can be readily 
ascertained. Tlie pressure Vv^hich causes the onset of yielding 
can be obtained using Tresca’s or von Mises' yield criteria 
( which are identical for the sphere) : 



where a is the yield stress in simple tension. This leads to the 
o ^ 

result that yielding (designated by an asterisk) begins when 



((iH-v) 
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a 

o 


( 7 ) 



Elas to-Plastic 

Here the material is elastic for r>c and plastic for r<c. In 
the plastic region the stresses are determined by Equation (6) - 
and the boundary conditions. The elastic strains in tiie plastic region 
are obtained from these stresses using Hooke's laws. It is 
assumed that strain hardening is absent, i.e., does not vary 
with the extent of plastic deformation. The stresses are continuous 
at the elastic-plastic boundary, r=c. 



The boundary conditions for /he elastic region are 


°0O' ''rr,= “o 


and 


0 = -P ..at r=b 

rr ext 


The solution for the elastic region is 


2 c ^ 1") ^ 

O ^ = -P , -4o f-3 (-3-1) 

rr ext 3 ob r 


and 


00 


2 r ^ h ^ 

» +4o 43 (t-3+1) 

ext 5 o b ^2r'* ^ 


The displacement in the elastic region is 

The general solution in the plastic region can be obtained 
from the equation of radial equilibrium, namely 


rdo 


rr 


dr-; ^ 


and Equation (6) . The result is 


a = 2 a il n r + D 

r r o 


where D is an integration constant which can be evaluated by the 
requirement of continuity of at r=c. Then we have 


-2o Jin— - P . - ~a (l-ra) 
ri o r , ext 3o^ b^^ 


( 9 ) 
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The local volume dilation in the plastic region is given by 


AV 

V 


(l-2v) 

B 



+a 


+ 0 


•)- (X-2v) 
rr-^ E 


(2a +3a 
^ o 



( 16 ) 


o 

The total change in volume of the plastic region (whose inner 
radius is constant a and outer radius is c) is therefore 


c 

(l^ _ 2v) (2a^ + 3a^^) 4fTr^ dr= 47 Tc^ u(c) . (17) 

a 

The value of is given by Equation (15) while u(c) 

follows from (13) with r set equal to c. It follows that 

(c/b)=Oo-|:>o • as) 


. = g-v) 

ext (1- 2v) 


(c/a) ^a_ 




For the fully plastic condition (labelled with double asterisk) 
c^b and (for v=l/3) 


* * 


P ^ = 2o K" - 2o ^nK, 
ext o o 


(19) 


and 


P ' - P . = 2a inK 

int ext o 


( 20 ) 


For the case where K is only slightly greater tlian one (the thin 

shell) we have from (19) P . = 2a 

^ K^l ext o 

as required since for a thin shell o^^ does not vary across the 
shell and the entire shell yields simultaneously as given by (8) . 



A ^ ^ 

The ratio vs K for v=l/3 is shown in Table 1. 


Table 1, Ratio of pressure needed to cause 
plastic flow throughout matrix over pressure 
needed for onset of plastic flow versus radius 


ratio . 

Based on Poisson 

ratio of one- third. 


K 

* A 

ext' 0 

P...MV2a 
ext o 

* A A 

P /P 

ext' ext 

1 

1 

1 

1 

2 

7.31 

0.71 

10.3 

3 

25.90 

0.68 

38.1 


As K becomes modestly large the pressure needed to cause plastic 
flow throughout the matrix is enormous relative to the pressure 
required for tlie conanencement of plastic flow. 

The extent of plastic flow for a given applied pressure and 
given K can readily be written in terms of f the fraction of 
the matrix which has undergone plastic deformation. Here 

I 

l“fj^ = (b^-c^)/(b^-a3) ^ (^21) 

Then (for v=l/5) from (18) 

^ )• (22) 

Curves of f^^^ vs ^"^ij 2 and 3 are shown in Figure 1. 

For the cases considered here, the strains everywhere are five 
percent or less (assuming 3(1- 2v) a^/E< 1/100) so the neglect of 
strain hardening is justified unless the matrix lias an unusually 
rapid rate of strain hardening. 



Alternate Method 


This problem can also be solved following the technique describe 
by Hill. The elastic stresses and displacements are identical 
as those derived above. In the plastic region, the equation of 
radial equilibrium again leads to tlie radial stress as in 
Equation (15). Now following Hill’s definitioji of tlie 
velocity o of a particle, we find 


du/dc 
1 - (9u/9f) 


(25) 


in the elastic region; and 


9o_^_2u 
9 r r 


- (l-2v) 


9c 


+^-^) (<^ ^^20 ) 
9 r^ ^ rr 0 0^ 


(24) 


in the plastic region. These velocities must be continuous at 
r=c,the plas tic- elastic interface. 

Using the equations for the stresses in the plastic region, 

(9) and (15), we can solve Equation (24) if we make the same 
approximations that Hill makes; i.e. neglecting second and higher 
orders of Oq/E, and neglecting 9u/9r witli respect to 1. Matching 
the velocities in the two regions gives the velocity in the plastic 
region, 

3(l-v)°o 2(l-2v)°o, c* c ®‘ext. , ,, 

B 7^ li b' 2^ 

O 


when is here written as a function of c. 

Since u is defined as dr/dc we can rewrite Equation (25) in 
terms of da/dc, tlic velocity of a particle at tlie inner boundary. 
With an incompressible material in the center da/dc is just zero, 
so we are left with an ordinary differential equation. 



_ {) _ 


Solvijig tills we find 


*’ext4 c>o(c/a)\F (26) 

where F is the constant o£ integration. We can calculate F 

A 

from the boundary condition that P =P for c=^a. This gives 

C JV L. C/ Jx C» 


For the fully plastic conditions 


A A 


ext 


U-v) 

TT^) 


a K^-2o £nK . 
o o 


(28) 


This result agrees witli Equation (19) in spite of the approximations 
used in obtaining (28) . 


CONCLUSIONS 

The pressure required to cause plastic deformation throughout 
a homogeneous elastically and plastically isotropic sphere surrounding 
a concentric rigid sphere is shown to be very large relative to 
the pressure required for the onset of plastic flow. In fact, tlie 
ratio of tliese pressures approaclies 1.5K^ for no strain hardening, 
where K is the radius ratio, for large K. Note that tlie ratio 
already approximates the large K limit for K=2 and K=3. For sucli 
values of K, strain hardening effects usually can be neglected, while 
for much larger K strain liardening would have to be taken into 



- :i 0 - 


account as would also nonlinear eJastic effects. 

For values of K of two or larger, tlie volume fraction of the 
matrix, which has undergone plastic deformation, increases nearly 
linearly with pressure above the critical pressure for the onset 
of yielding. 
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Figure Legends 

Figure 1. Volume fraction of matrix which is plastically 
deformed versus pressure for vaious K values. 

Here v=l/3. 

Figure 2. Ratio of internal pressure to external pressure vs 

applied pressure for K=2 and 3 for v=l/3. The single 
asterisk indicates the point of initiation of plastic 
flow. The double asterisk indicates the point at which 
• the sphere is fully plastic. The dashed line is the 
region in which Equation (20) holds. 






